Numerical simulation of small perturbation on an accretion disk 
due to the collision of a star with the disk near the black hole 

Orhan Donmez 

Nigde University Faculty of Art and Science, Physics Department, Nigde, Turkey 

(Dated: February 5, 2008) 

^O ' In this paper, perturbations of an accretion disk by a star orbiting around a black 

o, 

^^ . hole are studied. We report on a numerical experiment, which has been carried 

5h I out by using a parallel-machine code originally developed by Donmez (2004). An 

Oh, 
■^ ' initially steady state accretion disk near a non-rotating (Schwarzschild) black hole 

interacts with a " star" , modeled as an initially circular region of increased density. 

Part of the disk is affected by the interaction. In some cases, a gap develops and 

> 

Q> ' shock wave propagates through the disk. We follow the evolution for order of one 

^' 

CN ' dynamical period and we show how the non-axisymetric density perturbation further 

■^ , 

^5 , evolves and moves downwards where the material of the disk and the star become 

\o 

^^ ' eventually accreted onto the central body. 

^' 
Qh' When the star perturbs the steady state accretion disk, the disk around the black 



o 



hole is destroyed by the effect of perturbation. The perturbed accretion disk creates 



j^ [ a shock wave during the evolution and it loses angular momentum when the gas 

hits on the shock waves. Colliding gas with the shock wave is the one of the basic 

X" 

;_j ■ mechanism of emitting the X— rays in the accretion disk. The series of supernovae 

occurring in the inner disk could entirely destroy the disk in that region which leaves 
a more massive black hole behind, at the center of galaxies. 



I. INTRODUCTION 

It is widely accepted that the supermassive black holes (SBH) are located in the center of 
active galaxies and quasars (Begelman et al. 1984). The mass of SBH is usually estimated 
to be in the range M ^ 10^ — IO^Mq. Anomalous energy output of active galactic nuclei 
(AGN) may result from an accretion process. The accretion process involves matter which 
is attracted from the surroundings of AGN or comes from tidily disrupted stars passing 



too near by the SBH. The accretion disk is formed and eventually the matter falls into the 
black hole. The best evidence comes from studies of the central surface brightness of the 
nuclei, stellar velocity dispersion, spatial distribution of stars and chaotic time variations in 
X— ray emission (Young et al. 1978, Dressier et al. 1988, Dressier et al. 1990, Artymowicz 
et. al. 1993, Kato 2001). General relativistic effect is important consequences for the axial 
symmetry, stability and some of the phenomena of accretion disks. Artymowicz et. al. 
1993 estimates efficiency for rapid metallicity enrichment of the AGN and quasars based 
on star-disk interaction. The accretion disk around the massive black hole involves capture 
of low-mass stars from the host galaxy's nucleus by the accretion disk orbiting around the 
massive black hole. 

The dissipate interaction of star with accretion disk and main star was considered by 
different authors. The tidal interaction of star with central object, the massive black hole, 
was discussed by Hills (1975) and Artymowicz et. al. (1993). The interaction of star 
with an accretion disk was considered in the works by Vilkoviskij et al. (1981),Vilkoviskij 
(1983), Zentsova (1983), and Subr et. al. (2004). The gradual orbital decay of stellar 
trajectories is studied by Subr et. al. (2004). They have shown that how the stars become 
gradually flattened toward the accretion disk and identified parameters affecting the forming 
of structure during this process. The dynamical evolution of relativistic star cluster around 
a massive Kerr black hole with an accretion disk was examined by Ranch (1995). They have 
looked at the regime where the black hole dominates the potential and star-disk interaction. 
They have shown that there is an astrophysical plausible regime in which an accretion disk 
around the black hole in the AGN would dominate the dynamical evolution of star cluster 
through the star-disk interaction. Ranch (1995) and Artymowicz et. al. (1993) also declared 
that a series of supernovae occurring in the inner disk could entirely destroy the disk in that 
region which leaves a more massive black hole or neutron star behind. Ivanov et al. (1998) 
considered the hydrodynamical effects resulting from a collision among a supermassive black 
hole and an accretion disk orbiting another massive black hole which is bigger than the other. 
They studied these effects both analytically and numerically. They declared that the black 
hole-disk collision results in significant changes on disk structure, and gravity of the moving 
black hole induces a converging gas fiow in the disk. This converging gas helps to form a 
shock waves on the accretion disk. Subr et al. (1999) consider eccentric orbit intersecting 
disk once or twice per each revolution and they solve the equations for oscillating elements. 



They found the thermal radii of the orbits and time needed to bring the orbit into the disk 
plane as function of initial parameters, and they showed that previous simplified estimates, 
which is derived for the case of low eccentricity, remain valid within factor of two. The star 
disk interaction and gravitating disk around the black hole are reviewed within different 
approaches and with emphasis on the role of disk gravity by Karas et al. (2004). 

Theoretical prediction of star-disk interaction near a massive black hole is examined by 
Syer et al. (1991). Two different scenarios could happen when stars are orbiting on their 
eccentric orbit. The stars can pass through the accretion disk and create a hole or can be 
ground down by gravitational forces. The stars in galaxies could be ground down into a 
short period orbit during the star-disk interaction around the black hole. Any stars on a 
highly eccentric orbit, which are close to a black hole, would lose their potential energy. It 
may take a long time to be capture stars and grind down into an accretion disk. In Syer 
et al. (1991), they have considered mechanisms whereby stars can get on to orbits within 
Qcoii = GM/vl, where w* is the escape velocity from the star, M is the mass of the black hole 
and qcoii is major axis, without being destroyed. It might be thought that this could happen 
by a tidal capture process, if the stars get closer to the black hole than qt = (M/m^Y^^r^, 
they are tidily disrupted. A star captured into an orbit with pericenter at 3qt would undergo 
further tidal dissipations. The greatly increased geometric area and reduced binding energy 
suggests that such stars should be susceptible to large-scale stripping of mass on collision 
with the disk (Zurek et al. 1991). Within some critical radius Rent, collisions will then 
destroy the star and unbind the envelop mass within the red giant lifetime. The value of 
Rcrit will depend on the extent and properties of the disk as well as velocity and structure 
of star. If Rcrit and the central stellar density are sufficiently large, this mechanism could 
then provide an important source of fuel for the central black hole (Armitage et al. 1996). 

Different sorts of capture have also been suggested as formation mechanisms for binary 
stars. In these scenarios, some sort of process, which can dissipate kinetic energy, exist so 
that two unbound stars can become bound during a close encounter. These can be tidal 
capture where two stars through a close passage excite tidal modes in one another which 
absorb much of their relative kinetic energy and dynamic capture where three stars pass 
close together, exchanging energy so that one is ejected with a large part of the kinetic 
energy leaving the two other stars bound. The problems about these two models are that 
they require a high density of stars to make such close encounters statistically significant and 



thus probably can't form the majority of the overall binary population. It should be noted 
however that tidal capture is still the leading mechanism for formation of X-ray binaries in 
dense clusters 

In this paper, we numerically model the interaction of the accretion disk with star using 
general relativistic hydrodynamical code and assumed that the SBH forms such a binary 
system with a low mass star in an orbit inclined with respect to the plane of an accretion 
disk. The disk is thin and at the equatorial plane of the black hole but the model can be 
generalized to a more complicated model. We ignored the self gravity among the matter 
of accretion disk. The gravity between the black hole and accretion disk is set up by the 
Schwarzschild black hole. 

II. FORMULATION 

The General Relativistic Hydrodynamic (GRH) equations in Refs. Donat et al. (1998), 
Font et al. (2000), and Donmez (2004), written in the standard covariant form, consist of 
the local conservation laws of the stress-energy tensor T^'^ and the matter current density 



^^T>''' = 0, \7,J^ = Q. (1) 

Greek indices run from to 3, Latin indices from 1 to 3, and units in which the speed of 
light c = 1 are used. 

Defining the characteristic waves of the general relativistic hydrodynamical equations 
is not trivial with imperfect fluid stress-energy tensor. We neglect the viscosity and heat 
conduction effects. This defines the perfect fluid stress-energy tensor. We use this stress- 
energy tensor to derive the hydrodynamical equations. With this perfect fluid stress-energy 
tensor, we can solve some problems which are solved by the Newtonian hydrodynamics 
with viscosity, such as those involving angular momentum transport and shock waves on an 
accretion disk, etc. Entropy for perfect fluid is conserved along the fluid lines. The stress 
energy tensor for a perfect fluid is given as 

T^^^ = phu^'u'' + Pg'"'. (2) 



Here h is the specific enthalpy and 4- velocity m^ which may vary from event to event. It 
exhibits a density of mass p and isotropic pressure P in the rest frame of each fluid element 
The equation of state might have the functional form P = P{p, e). The perfect gas equation 
of state, 

P = (r - l)pe. (3) 

The conservation laws in the form given in Eq.(^ are not suitable for the use in advanced 
numerical schemes. To carry out numerical hydrodynamic evolutions such as those reported 
in Font et al. (2000), and Donmez (2004) and to use high resolution shock capturing schemes, 
the hydrodynamic equations after the 3+1 split must be written as a hyperbolic system of 
first order flux conservative equations. We write Eq. (^ in terms of coordinate derivatives, 
using the coordinates (x° = t,x^,x'^,x^). Eq.(IT)) is projected onto the basis {n^,{-£j)^}, 
where n^ is a unit timelike vector normal to a given hypersurface. After a straightforward 
calculation, we get (Donmez (2004)), 



dfU + diF' = S. 



(4) 



This basic step serves to identify the set of unknowns, the vector of conserved quantities 
U, and their corresponding fluxes FiU). With the equations in conservation form, almost 
every high resolution method devised to solve hyperbolic systems of conservation laws can 
be extended to GRH. 

The evolved state vector U consists of the conservative variables {D,Sj,t) which are 
conserved variables for density, momentum and energy respectively; in terms of the primitive 
variables {p,v\e), this becomes 
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Here 7 is the determinant of the 3- metric 7ij, Vj is the fluid 3- velocity, and W is the Lorentz 
factor. 



The flux vectors F* are, 
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The spatial components of the 4-velocity m* are related to the 3-velocity by the following 
formula: u* = W{v'^ — (3'^ /a), a and /5* are the lapse function and the shift vector of the 
spacetime respectively. The source vector S for Schwarzschild coordinate is given by Donmez 
(2004) 
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More detailed information related with numerical solution of GRH equations using high 



resolution shock capturing schemes in 3D are given in our paper Donmez (2004). This ref- 
erence explains detail description of numerical solutions of SRH equations, Adaptive-Mesh 
Refinement, high resolution shock capturing method used, and solution of GRH equation 
using Schwarzschild coordinate as a source term. The code is constructed for general space- 
time metric with lapse function and shift vector and is fully parallelized to make optimum 
use of supercomputers which is necessary to achieve the numerical modeling of real astro- 
physical problems, such as coalescing of the compact binaries and accretion disk around the 
compact objects. 

III. THE PHYSICAL MODEL AND CONSTRUCTING THE INITIAL 

ACCRETION DISK 



The structure and the dynamics of the accretion disk remain quite mysterious. The disk 
is not directly observable because the resolution of current telescopes is still insufficient. It 
is primarily studied at short wavelengths (UV, X and gamma rays). But short wavelength 
spectra give information only on the internal regions of the disk (scale of the micro-parsec) , 



very close to the black hole. One suspects that at these distances, the mass of the disk 
(generally regarded as small) starts to play a role on its own dynamics, and thus on its 
evolution and its structure. One then expects very particular effects, such as the generation 
of gravitational instabilities (spiral waves) 

The spiral structure on an accretion disk is inevitable in case of unstable situations. 
To simulate real accretion disk structure, the larger disk is needed. Despite large body of 
numerical work, there has been relatively little progress in the quantitative measurement 
of the effective transport of spiral shock in hydrodynamic simulations. This is due in large 
part to the fact that it is easy to set up a hydrodynamic simulation of a Keplerian accretion 
disk, but it is very hard to do it accurately enough to measure to relatively small transport 
expected due to spiral wave. To define structure of disk, one needs high numerical accuracy, 
a large range in radii, and long simulations covering several orbital periods. In order to 
watch the behavior of accretion disk and reduce the complexity, we restrict our attention to 
close to the black hole and put apply some restrictions explained below . 

A. Method and Assumptions 

To numerically model the accretion disk dynamics using the parallel-machine code, which 
can handle strong shock cases, we have made some assumptions. The computational model 
consists of two components: A central massive black hole and associated accretion disk lying 
in the equatorial plane of the hole with captured star. We presumed that the disk around 
the black hole is in the steady state condition. To restrict our attention to the influence of 
star-disk interactions on the evolution, and to keep the problem numerically tractable, the 
following assumptions on the calculations are made: i) the central black hole dominates the 
gravitational potential; ii) the accretion disk is thin and in steady state, which is created 
by numerical model; iii) the stars impact the disk supersonically; iv) captured star-disk 
interactions dominate the evolution of the stellar orbit. The thin accretion disk models 
describes an advection-dominated flow in which the radiative cooling is inefficient and most 
of the dissipated energy is advected into the black hole. In fact, the first constraint merely 
requires that the mass of the black hole be much larger then the mass of the cluster and disk. 
We neglect the self gravity. Outher boundary of numerical domain is chosen maximum in 
radial direction, Rmax = 30M, M is mass of the black hole, in where gravity of the central 
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massive black hole is too dominant compared with self gravity of disk. The validity of the 
assumption that star-disk collisions is dominating the evolution depend on both the density 
of the disk and the density of stars in the cluster. The problem of disks around the black 
hole is discussed within the general relativistic regime. Whereas general relativity is mainly 
relevant for the inner parts of disk (within lOOM Schwarzschild radius), Newtonian regime 
is adequate outside this radius. In the rest of the paper, distances in numerical calculations 
and time depend on the mass of the black hole M. 

B. Modeling a steady state initial accretion disk 

The numerical approximation of captured star by disk in the close binary system are 
modeled and evolved by solving the GRH equations and using fully GRH code (Ref. Donmez 
2004), run in two-dimensional spherical coordinates (r, 0) centered on the accretion disk, 
called equatorial plane. Theoretically, it is known that when the stars are orbiting on their 
eccentric orbit, they can be ground down or make a hole in the accretion disk around the 
black hole. The modeling of the perturbation needs an accretion disk, which could be in 
steady state or not, around the black hole. After the star is captured by disk, the region 
where star falls down during the evolution may not be affected by the star. On the other 
hand, steady state accretion disk can be used in the hole case which is created when stars 
are rotating in their eccentric orbit. Since there is no self gravity in our numerical code, the 
mass of the accretion disk does not effect the solutions. 

In order to simulate these problems, initial steady state accretion disk is computed at the 
equatorial plane (r — 0) of the black hole in the domains 3M < r < 30M and < < 2tt. 
The boundaries are periodic and the inner and outer radial boundaries are constructed to 
allow material to flow off the grid, called outflow, but not onto grid. This allows gas to fall 
into the black hole and out the computational domain when it reaches the inner boundary 
and the outer boundary in radial direction, respectively. An accretion disk is a structure 
formed by material falling into a gravitational source. Before applying any perturbation on 
the accretion disk, we start our analysis under the simplifying assumptions that the black 
hole has has steady state accretion disk around and it is numerically crated around the black 
hole choosing appropriate initial values. 

To create initial steady state accretion disk numerically, a sub-Keplerian velocity is chosen 
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as an initial velocity of gas because the gravitational forces will be balanced by the pressure 
and centrifugal forces on the accretion disk. Our simulations were run on uniform grid using 
256 points in radial and 128 in angular distances. For stability reasons, time steps change 
with changing sound velocity and velocity of matter. While initial disk simulation ran for 
over 11000 cycles, perturbation of accretion disk models ran for over 1510^ cycles. To create 
steady state accretion disk from the numerical simulation, the initial values of variables are; 
Pmax = 0.09, w^ = y/(M/r3)/^(l - (2M)/r) - l/r^, p = 0.1 and the radial velocity is zero 
outside the last stable circular orbit. The numerical simulation is performed until a steady 
state accretion disk created, seen in FigO The material in the quasi-steady state disk would 
keep on going in circular orbits and there would be no release of energy after the formation 
of the disk. Taking this as our initial condition, we then study the evolution of star-disk 
and gap-disk interaction close to the black hole using fully GRH code. 



IV. NUMERICAL RESULTS 

The dynamical evolution of captured relativistic star clusters around the Schwarzschild 
black hole with an accretion disk is examined in the regime where the black hole dominates 
the potential and star-disk interactions dominate the evolution. To take into account the 
interactions among the satellite and the disk matter due to general relativistic effects which 
are constituted by central object, numerical solution of the fully general relativistic hydrody- 
namical equations is needed. Theoretical prediction of star-disk interaction near a massive 
black hole is examined by Syer et al. (1991) and Subr et al. (1999). Two different scenarios 
can happen when the star is orbiting throughout its eccentric orbit. When the star can pass 
through in the accretion disk, it creates a gap and a hole in some region or it can be ground 
down by gravitational forces and in this manner, the accretion disk can be perturbed. 

Stars on highly eccentric orbits, which approach close to the black, would lose their 
orbital (potential) energy and angular momentum due to dynamic friction. The rate of 
star-disk captures is controlled by the rate of repopulating of orbits that intersect the disk. 
The captured stars at this scenario may give a perturbation and change dynamics of the 
accretion disk. Because there is considerable uncertainty in the actual value of the disk 
surface density and since the cluster consists of various stellar types ranging from compact 
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FIG. 1: Creation of the initial accretion disk during the time evolution. Top panel: density of 
the initial accretion disk vs. radial coordinate. The accretion disk reaches to quasi-steady state at 
t = 1100.5-/Vf. Bottom panel: radial velocity of gas during the creation of accretion disk. Radial 
velocity goes to zero outside the last stable orbit, around 6M, when disk reaches to steady state. 
This disk is used as an initial condition during the rest of the paper in our numerical calculations. 



stars to giants, there is also uncertainty in the star density that perturbs the accretion disk. 
Here, we have chosen a certain size of star, defined in section HVAI and it would be possible 
to have this size of star, and it can be entirely embedded in the disk when the inclination 
angle drops to critical value(Subr et. al. 2004). If the stellar mass is large enough to create 
a gap close to the black hole, the star becomes tidally coupled with the disk and the matter 
of the disk is sinked into the gap. Finally, these types of interaction perturb the accretion 
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disk and could create non-linear phenomena. 

This section describes two different GRH simulation varying the types of perturbation, 
grounded down star and gap on the steady state accretion disk using actual size of the black 
hole and a small perturbation which is represented by add density on the initial disk. 

A. Grounded Star 

To model the grounded star onto an accretion disk numerically, a perturbation is applied 
on the steady state accretion disk created in 1111 Bl During the interaction of star with accre- 
tion disk and central object, the star can be grounded down by forces which are produced 
by gravitational, centrifugal, pressure and tidal forces between captured star and the black 
hole. The grounded star can be represented as a high density region on the accretion disk 
and initial condition for the star is given as follow: ps = 0.1 located at 28M < r < 30M 
and < < 0.1. The density of star is %5 higher than the steady state density of initial 
disk and initial radial velocity, about expected velocity of star at the high density region, is 
given to star which is v^ = —0.01 (minus means star is falling to the black hole). The other 
variables are the same with initial steady state accretion disk. This is the scenario which is 
called a disk perturbed by a captured star 

Figs. El and 01 display the density at different times on the equatorial plane. The grounded 
star which has a higher density than accretion disk is moving toward the black hole under the 
general relativistic, centrifugal, pressure and hydro dynamical forces. The general relativistic 
future to the problem is given by the Schwarzschild black hole. While the perturbed gas is 
falling into the black hole on the equatorial plane, it creates one-armed spiral shock wave and 
shock wave propagates through the disk, shown in Fig. |2l The spiral wave moves through 
the disk causing the star and gas to clump up along the wave, called the density wave. In 
the early time of simulation, the shock wave created on accretion disk rotate in the same 
direction with the accretion disk and later, t = 1450. 3M, the spiral shock is settled down 
and starts to counter-rotate. The dynamic of disk during the later times of simulation is 
seen in Fig. (21 At the spiral wave front, where the density is higher, star formation may 
be enhanced in the vicinity of the black hole. Because the spiral wave travels through disk 
slower than the local Keplerian velocity, it carries negative angular momentum, and abrupt 
change or dissipation at the shock decreases the angular momentum of the orbiting gas on 
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disk. 

The shock wave created by the perturbation destroys the accretion disk. The spiral shock 
wave transports the angular momentum out of disk when the matter hits the spiral shock 
and it causes the gas falling into the black hole. The falling gas does not stop until reaching 
the black hole and the accretion disk is eventually destroyed by perturbation. Fig. |3] depicts 
mass of accretion disk during the numerical evolution. As can be seen from this figure, the 
mass increases after star is grounded down and then it exponentially decreases during the 
evolution. The ID cuts of density, orbital velocity and pressure vs. angular coordinate are 
shown at fixed radial distance r = 6M, Fig El The amplitude of shock created on accretion 
disk in the pressure and density plot is greater or equal to a factor of four. It is called the 
strong shock, converting the most of bulk kinetic energy to radiation and the thermal energy 
in the post shock material. Due to the shock waves, angular momentum and structure of 
disk is significantly changed. These parameters illustrate that moving shock wave is created 
on the accretion disk. The moving shock heats up material. We do not expect significant 
fraction of heated material which could have much effect on the overall dynamics of accretion 
disk. As a results of strong shock the matter loses its angular momentum and it begins to 
fall toward the black hole and a large amount of energy is released. 



B. Star Makes a Gap on an Accretion Disk 

As a second problem, the interaction of a gap with an accretion disk is modeled. When 
the star rotates in their eccentric orbit, it passes through disk and creates a gap in the 
accretion disk. It is theoretically explained by Syer et. al. (1991) and Subr et. al. (2004). 
In order to model this problem, the same initial steady state accretion disk is used with 
previous model, explained in llllBI 
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FIG. 2: The plots represent four subsequent moments of evolution of grounded star during the star 
disk interaction on the r — plane. The color corresponds the measure of density p. The captured 
star is located close to the outer boundary, seen at t = 43. 5M. The initial steady state disk is 
perturbed by star and disk matter moves in an unstable manner, and heats up accretion disk. As 
a consequence of this phenomenon, a one-armed spiral shock wave is created at later times. 

To model gap-disk interaction, gap is created on the initial steady state accretion disk, 
numerically, taken off the matter at 28M < r < 30M and < < 0.1. The removal of disk 
gas leads to the formation of a surface density depression that is continually replenished by 
the pressure forces and tidal interaction between gap and accretion disk. As consequence 
of tidal interaction, the spiral shock wave is created, seen in Fig. Fig. IHl illustrates 
the dynamical change of 2D representation of disk at different snapshot. Once the gap is 
created on disk closes to outer boundary, matter on the disk tries to fill out the gap and as a 
consequence of the filling, unstable manner is presented on the disk. This unstable manner 
forms one-armed spiral wave extending from center to the edge of the disk. Because of the 
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FIG. 3: Same as Fig|21but representation of density at later times. The spiral shock rotates in a 
counterclockwise direction while accretion disk is rotating in clockwise direction. 



spiral wave, surface density of disk exponentially decreases and hence, disk will not be able 
to survive. 

Fig. [TJillustrate the ID cut of density, pressure and orbital velocity vs. angular coordinate 
at r = 6M and t = 3740.87. The shock is numerically observed and rotating in the angular 
direction. The produced shock waves are similar and donate same physical behavior those to 
which are found from the grounded star simulation. These parameters illustrate that moving 
shock wave created on the accretion disk is called the strong shock. The shock waves first 
rotate in the same direction with accretion disk and then it settles down and counter-rotates 
to the disk. The moving shock heats up material. We do not expect significant fraction of 
heated material which could have much effect on the overall dynamics of accretion disk. As a 
results of strong shock, the accretion disk loses its angular momentum and the matter begins 
to fall toward the black hole undergoing a transition from subsonic to supersonic speed when 
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FIG. 4: The normalized mass of accretion disk during the time evolutions for grounded star case. 
Left panel shows the mass of disk during the hole simulation. The peak represents mass of disk 
in the moment of star captured. Spiral shock wave can act as a local sink for angular momentum 
and mass of disk decreases drastically, clearly seen in right panel (zoom of left panel around the 
small mass). This comes up with vanishing of accretion disk 

the disk material hits the shock waves. The results of perturbing disk by star have same 
behavior with making the gap on the accretion disk but the times for creating shock and 
locations are occurred at different times and different place, respectively. Creation of these 
phenomena depend on the gravitational, centrifugal, tidal and the hydrodynamical forces 
created on the disk around the black hole. 
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FIG. 5: The density, orbital velocity and pressure vs. angular direction for the grounded star 
model at r = 6M in ID. It is plotted at t = 2886M. Strong shock appears and it is moving around 
the black hole. 

V. CONCLUSION 



We have considered the general relativistic and hydrodynamical effect by solving general 
relativistic hydrodynamical equations resulting from perturbing accretion disk around the 
massive black hole. The theoretical predictions have shown that there is an astrophysical 
plausible regime in which an accretion disk around the massive black hole in an AGN or 
other galactic nucleus would dominate the dynamical evolution of a star cluster, through 
the process of star-disk interactions. We now discuss the general relativistic astrophysical 
consequences of star-disk interactions. As should be clear from the numerical simulation in 
this paper a lot of problems still remain in the context of binary formation. Many competing 
explanations for the most basic observations still exist and it is hard to distinguish between 
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FIG. 6: The plots represent four subsequent moments of evolution of perturbed disk by gap on the 
r — (j) plane. The color corresponds the counter of density p. The gap is located at close to outer 
boundary, seen at t = OM. Initial steady state disk is perturbed by gap and disk matter moves in 
an unstable manner, and heats up accretion disk. As consequence of this phenomenon, one armed 
spiral shock wave is created at the later times. The spiral wave moves at the same direction with 
disk at the early time of simulation. It rotates counter direction to disk after t = 2300M. 



them. But I think with the waist effort put into this field, theoretical predictions will be 
extended and better during the years, making it possible to distinguish really between the 
different theories and their outcome. This is however also needed as clearly, theories of 
binary formation are very important. Both since they probably must be applied to describe 
the actual formation of a majority of all stars but also because a lot of interesting phenomena 
like differences between star forming regions, bending of jets and formation of clusters can 
and must be discussed in connection with the formation and early evolution of binaries. 
The first numerical modeling results of perturbed disk are given in this paper. The star 
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FIG. 7: ID representation of Fig El for the density, orbital velocity and pressure vs. angular 
direction plotted at r = 6M and t = 2886M. The strong shock is presented and dramatic decrease 
or increase in the parameter causes to gas falling into the black hole. 

and the gap disk interaction close binary systems are analyzed theoretically and we made 
numerical modeling of a few of these predications after star opens a gap on an accretion disk 
or captured by an accretion disk. The numerical results show that regardless of whether stars 
is ground down or makes a gap, the star-disk interaction destroys the steady state disk around 
the Schwarzschild black hole. During this process all the matters at accretion disk fall into 
the black hole at the center and the black hole is getting supermassive. It is the important 
physical results. This may explain how the black holes are massive, 10^ — IO^Mq, at the 
center of galaxy and some of the mysteries in the accretion disk observations. Destroying 
an accretion disk and making massive compact objects (the black hole and the neutron 
star) is also declared by Ranch (1995). Although we star with steady state accretion disk 
to simulate this problem, non-steady state accretion disk case also gives same physical and 
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dynamical results. 

We stress that the above scenario seems to a remarkable degree model dependent, once 
density of star, size of captured star, initial angular velocity, location of perturbation in 
radial distance, and how often star can captured by accretion disk are specified. Our work 
may help to understand how the black hole at the AGN could be massive and why AGN is 
bright. 

In the star-disk and gap-disk interaction simulations, we have run the code on Cray T3E 
machine and Beowulf cluster paralyzing code using SCHMEM and MPI libraries. 
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